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Abstract 

Fixing the number of particles N, the quantum canonical ensemble imposes a con¬ 
straint on the occupation numbers of single-particle states. The constraint particularly 
hampers the systematic calculation of the partition function and any relevant thermo¬ 
dynamic expectation value for arbitrary N since, unlike the case of the grand-canonical 
ensemble, traces in the A-particle Hilbert space fail to factorize into simple traces over 
single-particle states. In this paper we introduce a projection operator that enables 
a constraint-free computation of the partition function and its derived quantities, at 
the price of an angular or contour integration. Being applicable to both bosonic and 
fermionic systems in arbitrary dimensions, the projection operator approach provides 
transparent integral representations for the partition function and the Helmholtz 
free energy F)v as well as for two- and four-point correlation functions. While appear¬ 
ing only as a secondary quantity in the present context, the chemical potential emerges 
as a by-product from the relation [Tat = F^+x - Fff, as illustrated for a two-dimensional 
fermion gas with N ranging between 2 and 500. 
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1. Introduction 


The calculation of the quantum mechanical partition function of N identical 
particles treated in the framework of the canonical ensemble remains a long-standing 
problem in many-body theory, even if the particles do not interact. The main difficulty 
hampering a systematic evaluation of Zat for moderate to large values of N originates 
from the particle number constraint that is to be invoked explicitly. In order to over¬ 
come this problem, we introduce a projection operator in section 2 which is capable 
of dealing with the particle number constraint for non-interacting particles (bosons, 
fermions) as well as systems of interacting particles complying with particle number 
conservation. However, the formal applicability to interacting particles is hardly useful 
in practice, because the eigenstates and the eigenvalues of the energy for such systems 
are rarely available. Although modern particle physics surely treats strongly inter¬ 
acting particles, it faces the necessity of applying approximations which, in essence, 
apply a variety of transformation techniques that reduce the problem to treating en¬ 
sembles of non-interacting particles. Thermal expectation values, based on statistical 
averages over ensembles of non-interacting particles, still provide the generic building 
blocks to set up perturbational and variational as well as other non-perturbative compu¬ 
tation schemes. Essential ingredients for such approaches are the partition function and 
the two- and four-point correlation functions characterizing systems of non-interacting 
particles. 

Keeping all this in mind, we believe it remains utterly relevant to consider a sys¬ 
tem of non-interacting particles and, therefore, we first examine its canonical partition 
function. As detailed in section 3, this results into a transparent integral representa¬ 
tion for the partition function as well as the corresponding Helmholtz free energy and, 
hence, the chemical potential of N non-interacting fermions or bosons. The integral 
representation also allows for a very simple derivation of a known recurrence relation 
for the partition function. 

For N harmonic oscillators in 1 dimension, the partition function could be ob¬ 
tained in closed form. The results are presented in section 4, for bosons as well as 
for fermions. 
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Section 5 contains a few numerical results related to the partition function and de¬ 
rived quantities of a finite size two-dimensional electron gas. Finally, the projection 
operator approach is applied once more in section 6 to derive generic expressions for 
the two- and four-point correlation functions. Some rather technical aspects are redi¬ 
rected to two appendices. 

2. The canonical partition function: a projection operator approach 

According to the nomenclature developed in the beginning of the 20th century, the 
statistical knowledge of a system in thermal equilibrium depends on the ensemble type: 
microcanonical, canonical or grand canonical. The canonical ensemble assumes that 
the exact number of particles in the system is known while its grand canonical counter¬ 
part merely requires that the average particle number be available. In theoretical studies 
of nuclear systems the number of particles is intrinsically dictated by the problem while 
for a great majority of solid-state systems only the average number of particles, in casu 
the density, is relevant. 

However, recent technological developments in nanoelectronics made it possible to 
control the number of carriers in nanometer-scaled devices, making the actual number 
of particles a more important parameter than the average number or density. Hence 
it would be desirable to export and extend theoretical methods developed in nuclear 
physics to various many-body formalisms commonly used to treat nanometer-scaled 
solids. A typical many body approach often starts with a short investigation of the 
non-interacting system, usually formulated in terms of creation and annihilation oper¬ 
ators. The use of these operators implicitly invokes a Fock space that, by construction, 
discards any reference to the number of particles whatsoever. However, if a descrip¬ 
tion with a fixed number is mandatory, one needs to introduce a projection technique 
that limits the Fock space to a subspace that corresponds to a fixed number of particles, 
while still allowing for a formulation in terms of the second quantization operators. The 
projection technique used for nuclear models can accomplish this task and is found to 
operate also for the second quantization description of a many-body Hamiltonian. Cor¬ 
respondingly, the number of particles is fixed and emerges as a fixed eigenvalue of 
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the number operator. After the projection one has to focus on the Fock subspace that 
is exclusively related to a fixed number of particles. In particular, the many-particle 
eigenfunctions of the projected Hamiltonian have to be calculated together with their 
energy spectrum and, afterwards, the probability density. 

Motivated by the above observations, we consider a fixed number N of of indis¬ 
tinguishable particles, fermions or bosons, described in the many-particle Fock space 
by a second-quantized Hamiltonian H. In order to preserve the number of particles, H 
has to commute with the particle number operator N. Consequently, due to [H, N] = 0, 
many-particle eigenstates of H can be found that simultaneously diagonalize H 


and N, i.e. 


H|(A/iv> = N|«Ai,iv> = = 0,1,2,3,... (1) 


Representing an arbitrary, allowable number of particles, the eigenvalues A of N are 
used to label the eigenstates as well as the corresponding £j;v. The index j co¬ 
vers all remaining, internal quantum numbers that are labeling Ej f^ for a fixed value 
of N. For the sake of notational simplicity, we have omitted below any dependence 
on spin components which, however, can be incorporated into the formalism when¬ 
ever required. Because H operates in Fock space without any a priory reference to the 
number of particles, thermodynamics is usually expressed in the grand canonical en¬ 
semble (GCE). Within this framework, the chemical potential emerges as a Lagrange 
multiplier regulating the average number of particles, rather than imposing a sharply 
defined value of N, as required in the canonical ensemble (CE). In order to overcome 
this problem, we propose a projection operator that extracts a A-particle Hamiltonian 
Hat out of H, while automatically invoking the canonical constraint of N particles. 

Let {ji/'/i.M)} denote the complete set of eigenstates with an integer, nonnegative 
eigenvalue M of the number operator N: 



( 2 ) 


and consider the operator 



( 3 ) 
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( 4 ) 


with the obvious properties P|, = P^f and P^ = P^^. One furthermore observes that 

1 

Pv \fn,M) ^ ^ J d0 |^«,m) = ^JV,M \^n,N) 

and hence 

Pv ^ \^n,M) - ^n,N |t^(7,v) ■ (5) 

n,M 

Therefore, P^f is a real projection operator which yields an eigenstate of the A^-particle 
subspace if it acts on an arbitrary state of the entire many-particle Hilbert space. Con¬ 
sequently 

Hv = PvH Pn (6) 

is the A^-particle Hamiltonian, extracted from the Hamiltonian H in the many-particle 
Hilbert space. Note that the position representation of Hat in principle coincides with 
the N-particle Hamiltonian of first quantization, as can be inferred from the algebraic 
treatment given in Ref. [ 1 ]. Although similar projection operators have been introduced 
before in statistical physics [2], nuclear and high-energy physics [3, 4, 5], we are not 
aware of its practical use as a particle number regulator in quantum statistics. 

It is tempting to immediately suppose that the partition function for thermodynam¬ 
ical equilibrium is given by 

ZwOS) = Tr(e-^^«) (7) 

with some typical derived quantities as the Helmholtz free energy Tjv(j8) and the internal 
energy with yS = lUksT), where is the Boltzmann constant (ks = 1.3806568 x 

and T the temperature in Kelvin: 

FN(fi) ^-UniZj.m ( 8 ) 

W) = -4ln(Z^0S)). (9) 

Although correct, these equation should be handled with care. Thermal equilibrium 
means that the internal energy Un is stable in time, and p (and hence T) is in essence 
a Lagrange multiplier for imposing that condition, rather than a given quantity. The 
internal energy Un is the fixed quantity. Because of the technicality of this question, 
the correct interpretation of the principle of maximum entropy [6, 7, 8] in thermal 
equilibrium is treated in Appendix A. 
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This approach is clearly consistent with (A. 12) under the condition (A.6). Defining 
the generating function 

GO6,0) = Tr(e-^^e‘'^®), (10) 

and comparing it with the definition (3) of FV, we obtain as the A-th coefficient 
of the Fourier series that represents GijS, 8), 

G(p,0)s-'^^d0. ( 11 ) 

3. Partition function of non-interacting indistinguishable particles 

As already argued in the Introduction, the general formulation of the previous sec¬ 
tion, though valid also for interacting particles, is of limited practical use. Quantum 
statistics of non-interacting particles on the other hand, still provides the basic ingredi¬ 
ents for most approximative treatments of interacting particles. Therefore, we first con¬ 
centrate on the partition function of non-interacting particles with supposedly known 
eigenstates and energy levels. The Hamiltonian H and the number operator N can then 
be expressed in terms of the single-particle energy spectrum e^., where k denotes any 
set of generic quantum numbers properly labeling the single-particle energies and the 
corresponding eigenfunctions: 

H = N = = (12) 

k k 

where the creation and destruction operators cl and satisfy appropriate (anti)commu- 
tation relations, i.e. 

[cl, c[] = 0 = [c^, c^,], [cj,, c[] = 6k,k’ for bosons, (13) 

{“it ’ <^1'} = 0 = {c^, c^,}, ^k,k' for fermions. (14) 

This means that any particular energy in (1) takes the form 
Ei,N^Yj'ikh with 

k k 

the integer occupation numbers being restricted to 0 and 1 for fermions while rang¬ 
ing between 0 and infinity for bosons. Keeping the total number of particles fixed is 
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prohibitive [9] for writing as Ha exp(-y6«^e^,). As can be found in many text¬ 
books, e.g., in Ref. [10], the standard approach to remedy this problem involves the 
construction of all cyclic decompositions of the particle permutations, which turns out 
to be a tedious task. Use of the projection operator greatly simplifies this conditional 
summation. Elementary operator algebra enables one to work out (10) explicitly, yield¬ 
ing 

G({5,0) = Tr(nA = Ha and thus 

nz exp((i6'-/3e^)n^). (16) 


Summing from 0 to oo for bosons, and from 0 to 1 for fermions, readily gives 

G(y6,0) = j~~[ (l - exp(i0forbosons, (17) 

A 

G(Ji, ^) - j~~[ (1 + (i® “ ^^a)) fermions, (18) 

A 

which (less transparent but more compact) can be abbreviated as 


G(P, 0) = j~[ (1 - ^ exp (i0 - ySej.)) ^ with 

A 


^ = +1 


for bosons, 
for fermions. 


(19) 


Filling this out in (11), it should be noted that the angular integral can equivalently be 
expressed as a complex contour integral along the unit circle 


= = ( 20 ) 

The generating function G(J5,z) is analytic everywhere for fermions (^ = -1) whereas, 
for bosons, the region | z | < 1 would merely contain an isolated singularity at z = 1 if 
the single-particle ground-state energy were vanishingly small. In order not to intro¬ 
duce spurious poles, all boson single-particle eigenenergies should be strictly positive. 
This can always be realized by an energy shift resulting from a gauge transformation. 
This ensures that, inside the unit circle, the integrand of Eq. (20) has a single pole at 
z = 0, whence 


Zn(/3) = ^ lim 
A^! 


d^G(J5,z) 

dz^ 


( 21 ) 
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Obtaining first the derivative of In G(J5, z) to get 



( 22 ) 


we apply Leibniz’ rule to take the {N - l)-th derivative of Eq (22) for ^ 1 to arrive 
at 


Z„(je)=l, Zi08) = ^exp(-y6e,), 


k 





This recurrence relation is not new [11] [12], but its derivation from the contour in¬ 
tegral (20) is substantially simpler than what follows from a tedious analysis of the 
permutation group. It also enhances the confidence in the correctness of the projection 
operator approach. 

Given the occurrence of the variable z as a prefactor of the exponentials in 
Giji, z), it might be tempting to interpret z as a complex fugacity in analogy with the 
real fugacity expfySqccE) appearing similarly in the grand-canonical partition function 
and the Bose-Einstein and Eermi-Dirac distribution functions. However, a safer inter¬ 
pretation could lie in the comparison of the CE and the GCE: whereas the latter sets 
the chemical potential to fix the average number of particles arising consequently as a 
weighted sum over all particle numbers N, the CE fixes N and is therefore bound to 
integrate over all relevant “complex fugacities”. 

To clarify this point, we extend the unit circle in (20) to another circular contour 
with radius u > 0: 



(24) 


The fact that this expression is independent of u implies dZf^(P)ldu - 0. Because 
of Eq. (22), this means that 
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d6» = 0. 


(25) 








The above sum rule for the CE can not be satisfied by mgce, the value of m = exp(yS[^i) 
that solves the transcendental equation for the GCE, i.e. 


k 

Consequently, in the light of the CE, Eq. (26) should be considered an approximative 
equation, usually obtained from a saddle point method. The latter amounts to maxi¬ 
mizing the factor G(J5, u)/u^ in the second line of Eq. (24), where it is expected that 
Zffiji) ~ G(J3 ,ucce)Iu^ce ^ good estimate of the free energy. And indeed, the 
Helmholtz free energy then becomes 


FN(fi)^ 


^ ^ In ^1 - e^hi ft) j + ]\f^i for bosons. 


i ^ In ^1 H- e^^^‘ j + N\i for fermions. 


(27) 


Eor sufficiently large N this is consistent with the familiar assumption p 
Fn(J3), since then Pat+i ~ bv- But the present derivation clearly shows how and why the 
standard transition from the CGE to the CE is an approximation. A correct treatment 
of the CE has to deal with the angular integral or, equivalently, the complex contour 
integral for (or its equivalent recurrence relations). 


4. Indistinguishable harmonic oscillators in ID 

Until now, closed form solutions involving indistinguishable particles are barely 
available, even if they are not interacting. As an exception, however, we illustrate the 
case of non-interacting bosons and fermions collectively moving in a ID harmonic 
potential and sharing the well-known single-particle energy spectrum 

= k = 0, l,...,oo. (28) 

The countour integral representation (20) or, equivalently, the derivative rule (21) re¬ 
lates to the generating function G(J3,z). In the present case, the latter is given 
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by 


/ oo 



1 _ 7 p-yS^a>(^+l/2) 


k=0 


=n 


for bosons. 


for fermions. 


(29) 


Direct evaluation of the A^-th derivative of G(J3,z) seems quite a formidable task, if 
possible at all. However, two mathematical identities derived by Leonhard Euler and 
nowadays emerging as corrolaries of the ^-binomial theorem [13, 14] (see also Ap¬ 
pendix B) are found to solve the problem. According to the identity (B.l), the infinite 
product for bosons can be written as a convergent series for |z| < 1. Similarly, the 
identity (B.2) can be used for the fermionic case. The result is 


G0S,z) = 


-phcjk 


CO n ^ 

n=l k—l 

oo n ^ 

1 - 1 - Y\ __ 


-phcjk 


for bosons. 


for fermions. 


(30) 


In accordance with Eq. (21) the coefficient of in the above series is the partition 
function for N oscillators: 


ZwOS) = 


-Nphwll 


-fifiLijk 


^-N^I3tioj/2 






k=0 


for bosons. 


for fermions. 


(31) 


Having determined the partition function, one may easily find the Helmholtz free en¬ 
ergy 


Fn(/3) 


1 ^ 

iyin(l 

pit ^ 


■ H- 


1 

-NHlo 

2 

1 9 

-N^fio) 


for bosons, 
for fermions. 


(32) 


Complying with the standard definition Pat = of the chemical potential, one 

thus readily obtains 


ln(l - Q-iN+k)l}tioj\ 
fiAf(^) = - — - + 


ho) 


A + - I fiw 


for bosons, 
for fermions. 


(33) 
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clearly depending on both N and f}. Only for sufficiently large N, more precisely for 
^-Nphw ^ j-jjg logarithmic term can be neglected, such that 






Tiu 


N + 


For the internal energy, one finds 


t/ivOS) = ^ 

k=\ 


kfioj 


Q^hcLjk _ 2 


for bosons, 

Ha> for fermions. 

-NHo) for bosons. 


(34) 


1 


-N^Tioj for fermions. 
2 


(35) 


As discussed in Appendix A, this equation should be considered as a transcen¬ 
dental equation, determining p for given L)v- Nevertheless, it is common practice to 
express Un as a function of (3, in which case it would just take a simple rotation of the 
corresponding UNifi) curve to obtain the requested /3(Un) relation. 

However, it is more instructive to look at the specific heat dL^v/dr, usingy6 = 1/ (ksT): 


dU, 


— = /tfi y 

dr ^ 


A=1 


(k/rf e^’/^ 
(e^/^ - 1)^ 


with 


ksT 
fiu ’ 


(36) 


which holds for both fermions and bosons. It is clear that dL)v/dr —» Nk^ as ex- 

T —»oo 

pected. The important point to be emphasized is that the relation between Un and the 
temperature is in general not linear. The convergence to the classical limit Un Nk^ T 
even slows down with increasing N. 


5. A two-dimensional electron gas 

The formal expressions for the partition function obtained in sections 3 clearly 
show that any practical investigation of statistical physics within the framework of 
the canonical ensemble is bound to deal with the angular integral or, equivalently, the 
complex contour integral. It occurs in the expression of all thermodynamical quanti¬ 
ties (partition function, free energy, specific heat, ...), either by direct evaluation or 
by conversion into the equivalent recurrence relations. Analytical results can only be 
expected for an extremely small number of systems. The previous section gave such 
an example, but in general one has to rely on numerical methods. 


11 










Figure 1: Helmholtz free energy for Lx = Ly = 100 nm. 


As an illustration we quote the calculation of the free energy and the chemical 
potential for a free electron gas residing in a finite, two-dimensional rectangular area 
0 < X < 0 < y < Ly, while imposing periodic boundary conditions on the single¬ 

electron wave functions and taking the single-electron energy spectrum to be 


ti^ 

2m 


2nn, 


2mi^ 

~ 


tijc, riy — 0,±l, ±2,... 


(37) 


with OT the electron mass. 

In view of possible practical applications, for example for the electron gas in the 
inversion layer of a MOS field-effect transistor, we have fixed and Ly to be 100 
nm, whereas the ambient temperature is assumed to be 300 K. Extracted from the 
recurrence relation (23), the Helmholtz free energy is plotted in Fig. (1) versus the 
number of electrons N, up to A = 520 beyond which sign changes (due to ^ = -1 for 
fermions) made the recurrence relation unstable. It turns out that, for fixed Ly = L,., the 
free energy attains a minimum for a particular value of A - in the present case around 
A = 374 - which corresponds to = 0 or, equivalently, the absence of energy cost 
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2DEG chemical potential 



Figure 2: Canonical and grand-canonical chemical potentials for Lx = Ly = 100 nm up to = 520. 


when a single particle is to be added or removed. On the other hand, a typical value 
of the areal electron concentration in a MOS capacitor operating at room temperature 
is 10*^ cm“^ which, for - Ly = 100 nm, corresponds lo N ^ 100 and, hence, to a 
negative chemical potential. 

Unlike the grand-canonical chemical potential pccE that corresponds to the 
thermodynamic limit N —» oo,L^, L,, ^ oo, whilst «s = NfiL^Ly) remains finite, can 
be calculated analytically from 


Pgce - 


-in (exp I 





(38) 


For the sake of comparison we have plotted both and pccE versus N in Fig. 2 and in 
Fig. 3. For relatively large values of N, say N > 30, Pai and pocE are about equal. On 
the other hand, the expression (38) for pccE is only valid in the thermodynamic limit. 
It fails to characterize electron ensembles with < 30 which can, however, be handled 
by the canonical formalism yielding pai. 

Moreover, in the case of more complicated fermionic systems, such as the 3DEG, 
the transcendental equation (26) expressing the (average) number of fermions in terms 
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2DEG chemical potential 



Figure 3: Canonical and grand-canonical chemical potentials for Lx = Ly =100 nm up to = 30. 

of Pgce generally can no longer be inverted analytically, while the computational 
scheme that yields remains unaltered. 

6. Correlation functions 

So far we have concentrated on the projection operator approach for obtaining the 
partition function and derived quantities of indistinguishable particles, with particular 
attention to non-interacting particles. Two specific examples were worked out. But, 
as already mentioned in the Introduction, not just the partition function and its derived 
quantities but also single- and two-particle correlation functions of non-interacting par¬ 
ticles are instrumental to perturbative and variational methods that are commonly en¬ 
tering approximative treatments of interacting particles. In the present section, it is 
shown that the projection operator is also well equipped to calculate these quantities. 

Single- and two-particle correlation functions - also referred to as two-point and 
four-point functions - typically provide a signature of the correlation between particles 
that are spatially separated. For the sake of notational simplicity, positions in space are 
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denoted by x and x' which, however, should not at all be regarded as a limitation to 
strictly one-dimensional systems. 

Considered quantum statistical averages, correlation functions are conveniently ex¬ 
pressed in terms of field operators ili{x) satisfying typical (anti)commutation relations 
{i/'(x), ij/^(x')} - 6{x - x') etc. 


6.1. Single-particle correlation functions 

Adopting once again the canonical ensemble framework, we can reinvoke the above 
defined projection operator FV to calculate the single-particle correlation function (two- 
point function) for an ensemble of N particles from 

S 2 (x, x') = (lA^(x) Iff(x')) = ^Tt(pN exp(-je H) f^x) fix')). (39) 

Z.N 

In most cases of interest, both the single-particle energies and the corresponding 
single-particle wave functions (pk{x) are supposed to be explicitly known, the latter 
constituting a complete, orthonormal basis. Hence, it proves convenient to expand the 
field operators as 


k' k 


the creation and destruction operators thereby appearing as expansion coefficients. 
Substitution into (39) yields 


S 2 (x, x') = ^ (pKx) (pk'ix') (41) 

kk' 


with 


= ^Tr(P^exp(-/lH)c/c,,) 

T 

d0 exp(-iA0)Tr(exp(i0N 

T 

-i— r d0exp(-iA0)Att,(0). (42) 

ZttZn J-^ 


In order to evaluate the trace 



A«'(0) = Tr(exp(i6>N -ySH)c/c^,), 


(43) 
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we first exploit its invariance under cyclic permutations to get 

^kk'(0) = Tr (c^, exp (i0 N - H) c/). (44) 

At this point, we insert the identity operator exp(i6( N -/i H) exp(-i6( N +y6 H) under the 
trace, in front of c^,, 

^kk'(0) - Tr(exp(i0N -;6H)exp(-i6N +y6H)cj,, exp(i0N (45) 

and apply the operator identity 

exp (-i0 N + jS H) c^, exp (i0 N - /i H) = exp (i0 - c^, (46) 

that proves valid for non-interacting particles. Indeed, given the Hamiltonian H = 
Yjk Cf. both exponents in exp [-\0 N + y6 H) c^, exp (i0 N - >6 H) are found to factor¬ 
ize while c^, commutes with each factor but the k'-i\\ one. Consequently, all factors 
appearing in the right exponent other than the ^'-th one can be shifted to the left so as 
to neutralize their inverse counterparts. Hence, we are left with 

exp(-i0N +/3 H)c^, exp(i0N -/3H) = exp(zc^,c^,)cj,, exp(-zc^,c^,) , (47) 

where z = - 10. Differentiation of M/t'(z) = exp(zc|,c^,)cj,, exp(-zc|,Cj,,) with 

respect to z yields a first-order linear differential equation 

- exp [zclc,) [4c,„ c,,] exp {-zclc,) 

= - exp exp (-z^Cj.,) = -ut{z) (48) 

to be solved with the boundary condition Uk>{Q) - Q'. The trivial solution Uk>{z) - 
exp (-z) Ck’ immediately leads to the operator identity quoted in Eq. (46). 

As a result, we obtain 

Kkk’(0) = exp(i0-/ie^,)Tr(exp(i6'N -/?H)c^,c/.). (49) 

Exploiting Cj^,cl - we can rewrite the above results as 

Kkk'{0) = exp (ie - ySe^„) \6kk'G(J5,0) + ffr (exp (i0N -y8 H) c/c^,)] 

= exp (ie -ySe,,) [dkk'Gifi, 0) + ^Att-(0)]. (50) 
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Hence, the trivial solution reads 


Att'(0) - 6k’k- 


G(J5,0) 


exp(/3e^-i6»)-^ 

In turn, the expression for the single-particle correlation function simplifies to 


(51) 


S 2 (x,x') = 'y 4>l(x)(pk(x') f d0 exp(-iN0) - - ^^2) 

IttZjv Y exp (ySe^ - i 0) - ^ 

with the particle density n(x) = S 2 (x, x) emerging as a particular case. 

6.2. Pair correlation functions 

Introducing the pair correlation function (four point function) as 

54(x, x') = (fHx) fix) f\x') fix ')), (53) 

we first write S 4 ix, x') as 

S4ix, x') = if^ (x) fix')) 6ix — x') + if^ (x) f\x') fix') fix)) 

- nix) Six - x') + if^ix) f\x') fix') fix)). (54) 

Expanding again all field operators in the complete set [ftix)), we obtain 

{f\x) f\x') fix') fix)) = EE ^*qix) f‘g,ix')(l)k'ix') fkix) {cl cl, Ci^, cf) . (55) 

kk' qq' 

A lengthy but straightforward calculation involving another application of the operator 
identity (46) and the commutation relation [c^, cj c^,] = dkqcl, + fSkq'cl leads to 


(cl c^, c,,c,) — - 

\ g q' k' k' 2nZk 


■ i^^kqSk’q' + Skq'Sk'q )£- exp i—iN0) 

Gif,6) 


(exp (j3e^ - 10 ) - f) (exp (/3e^, - w) - f) 
Correspondingly, the pair correlation function is given by 


■ (56) 


{f\x) f\x') fix') fix)) = y [f\(j)kix)?'\fk'ix')f- + (pl,ix) flix')(pk'ix') fkix)) 

Gif},0) 


kk' 


f 


X I d6exp(-iA0) 


(exp (/?e^ - 10 ) - f) (exp {fe^, - w) - f) 


(57) 
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Appendix A. Principle of maximnm entropy 


Consider the density operator 

pN ^ 'y^Pi,NWj,N){>l/j,N\ with (A.l) 

j j 

where is the probability that state of the A-particle subspace is occupied in 
thermal equilibrium, i.e. with a fixed ensemble average Ut^ for the energy [6, 7, 8]: 


Un - Tr^HwpA') - PinE^n- 


(A.2) 


Maximizing the entropy 


imposes 


Sn - -^fiTrlpArlnpiv) - -^b '^Pj.N^np 


'j,N 



^ Pf,N (in p/,iv + O' + y6£;',iv) = 0, 

/ 


(A.3) 


(A.4) 


where a and p are Lagrange multipliers for the normalization and the energy condition, 
respectively. Hence p^N - withZA^ = But, keeping in mind that Zat 

and p are in fact functions of the fixed value Un, a more careful notation is introduced: 


Pj,N{UN) - 


_ Q-pi )Ei,N 

'Z.n{Un) 


and therefore 


j 


(A.5) 


Un = 


1 

Zn{Un) 




Sn(Un - ks InZNiUN) + k^p{UN)UN- 


(A.6) 

(A.7) 


where the equation for Un is a transcendental equation which determines the Lagrange 
multiplier/?, and hence the temperature T if defined as /? = ipk^T). The Helmholtz 
free energy 

Fn(Un) = Un - , Sn(Un) (A.8) 

ksPiUN) 

then becomes, as expected: 


Fn(Un) =--^^\nZN(UN). 
PWN) 


(A.9) 
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At first glance all these result are familiar. Less familiar is a relationship between the 
entropy and the energy dependence of p. Differentiating with respect to t/^v, one 
obtains IAUn - IdU, and hence 



Using (A.6-A.8), this expression simplifies into 



(A. 11) 


showing how the entropy increases with increasing internal energy. In the classical 
limit, with p{U) - CPk^U) where C is the specific heat, the right hand side becomes 
C\n{UNIU q) - Cln(7'A^/7’o), consistent with the equipartition theorem. 

So far, it was shown that the projection operator approach is consistent with the 
standard interrelations between the thermodynamic quantities, all derivable from the 
partition function Zn and the (given) internal energy Un- No attention was paid to the 
relevance of the projection operator for the actual calculation of Zat, which becomes 
now the main topic of interest. Since = FV, and Z^ can be rewritten as Z^ - 

Tr(e-/^n«), 

one obtains with little effort from (A.5) 


Z^(t4)=Tr(e-^"*>Hp^), 


(A. 12) 


regardless whether the particles are interacting or not. Without the projection operator, 
this would be the grand canonical partition function, for which the chemical potential 
is required as a Lagrange multiplier to impose the average number of particles. The 
present approach is bound to work in the canonical ensemble with exactly N particles. 

Until this point, a purist notation was followed, emphasizing that thermal equilib¬ 
rium means that the internal energy Un is fixed, and that a Lagrange multiplier P{Un) 
is introduced in (A.6) to fulfill this requirement [6, 7, 8]. For practical purposes, this 
formal treatment is less appropriate. It is much easier to consideryS as a function argu¬ 
ment 


Z^(/3) = Tr(e-^^P^), 


(A. 13) 


which at the end of the calculations is connected to the internal energy via 


mP) = :^Tr(He-^^PJ = -Ain(z^ 08 )). 
Zn{P) '' ' dp 


(A. 14) 
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Appendix B. Two Euler identities 


Given two complex numbers z and q, with \q\ < 1, Leonhard Euler in the 18th 
century derived (amongst a variety of other mathematical insights) the following two 
identities: 


OO 


n 

n =0 


1 

1 - zq" 




1 

1 -q!^ 


(B.l) 


uo yo n ^ 

f[ (1 - zq") = 1 + Xi n 

«=0 «=1 k=\ ^ ^ 

In contemporary literature [13, 14, 15, 16], they are usually obtained as a by-product of 
more general theorems on ^-products and ^-series, which hinders a transparent deriva¬ 
tion. Therefore we propose an easily accessible proof, inspired by a strategy of Berndt 
[16]. Given a set of complex numbers a,b,q,z with |z| < 1, \az\ < 1, \bz\ < 1, \q\ < 1, 


define a function 


and calculate f(qz). 


fiqz) = Y\ 

n =0 


1 - 

1 - bzq"^^ 


CO 


/ft)=n 

n -0 


1 - azq’^ 
1 - bzq’' 


1 - azq" _ 1 - fe 1 - azq" 
^ I 1 - bzq’' 1 - az ^ i 1 “ bzq" 

«=1 n =0 ^ 


The latter can conveniently be rewritten as 


(B.3) 


1^/ft). (B.4) 
\ - az 


f{z) - f(qz) = bzfiz) - azfiqz). (B.5) 


Since /(z) is analytic wherever |z| < 1, we may assign a power series to it: 

OO 

fiz)^Yu^„z'', (B.6) 

n =0 

where Co = /(O) = 1 holds by construction of /(z). 

Substituting (B.6) into (B.5), we obtain 

OO OO 

2 C„ (1 - q") z” = ^ C„ (17 - aq") z"+'. (B.7) 

«=0 /i =0 

Clearly, the n = 0 term in the left-hand side of Eq. (B.7) vanishes, while its right-hand 
side may be rephrased by shifting the summation index 7z -i- 1 —> «, yielding 

OO OO 

2 C„ (1 - q") z" = ^ C„-i [b - aq”-^) z«. (B.8) 

«=1 n =\ 
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Identification of the coefficients of z" then leads to the recurrence relation 



(B.9) 


which can be solved with the help of Co = 1 to get 



(B.IO) 


Filling C,, out in (B.6) gives a generalization of the well-known ^-binomial theorem 



(B.ll) 


The Euler identities (B.l) and (B.2) emerge as special cases of (B.ll), corresponding 
respectively to the cases a - 0,b - \ and a - l,b - 0. 
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